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Abstract 

This  lecture  presents  the  role  of  Polarimetry  in  SAR  Interferometry.  A  general  formulation  for  vector  wave 
interferometry  is  presented  that  includes  conventional  scalar  interferometry  presented  in  the  respective 
former  lecture  as  a  special  case.  Based  on  this  formulation,  the  coherence  optimization  problem  can  be  solved 
to  obtain  the  optimum  scattering  mechanisms  that  lead  to  the  best  phase  estimates.  Comparison  with 
conventional  single-polarization  estimates  illustrates  the  significant  processing  gains  that  are  possible  if 
there  is  access  to  full  polarimetric  interferometric  data.  A  comparison  with  conventional  single-polarization 
presented  in  former  lectures  illustrates  the  significant  processing  gains  that  are  possible  if  access  to  full 
polarimetric  interferometric  data  is  possible.  The  strong  polarization  dependence  of  the  coherence  will 
be  addressed  and  the  analytical  solution  for  optimum  polarization  states  that  maximize  the 
interferometric  coherence  will  be  derived  and  applied  to  experimental  data.  These  improved 
interferogrammes  allow  an  improvement  of  the  accuracy  of  derived  DEM  products. 

The  introduction  of  a  new  coherent  decomposition  theorem  for  interferometric  applications  based  on  the 
Singular  value  spectrum  of  a  3  x  3  complex  matrix  allows  the  decomposition  of  polarimetric 
interferometric  problems  into  a  set  of  coherent  scattering  mechanisms.  As  a  consequence,  it  is  possible  to 
generate  interferograms  related  to  certain  independent  scattering  mechanisms  and  extract  the  height 
differences  between  them.  The  limitation  of  this  technique  is  the  existence  of  independent  scattering 
mechanisms  located  at  different  height  positions. 

To  explain  the  physical  origin  of  these  mechanisms,  a  coherent  electromagnetic  scattering  model  will  be 
established  which,  additionally,  can  be  used  to  establish  the  suitability  of  the  decomposition  algorithm 
for  solving  the  problem  of  estimating  the  location  of  the  effective  scattering  center,  which  is  a  critical  point 
in  the  physical  interpretation  of  interferograms.  However,  these  introduction  of  Polarimetry  in 
interferometric  processing  requires  that  fully  coherent  polarimetric  data  must  be  collected  in  order  to 
separate  the  scattering  mechanisms..  In  this  connection  polarimetric  Differential  Interferometry  will  be 
considered  also. 

The  phase  difference  between  the  optimum  interferograms  obtained  by  the  application  of  the 
algorithm  on  the  SIR-C  data  turned  out  to  be  strongly  correlated  with  the  actual  forest  height.  This  was 
a  major  result  indicating  the  potential  of  the  coherent  combination  of  polarimetry  and  interferometry. 

Introduction 

Polarimetric  SAR  Interferometry  was  a  first  step  in  the  abatement  of  the  scattering  ambiguity  problem  in 
the  height  direction.  By  combining  interferometric  and  polarimetric  techniques,  it  enables  the  separation 
of  different  scattering  mechanisms  within  a  resolution  cell  and  at  the  same  time,  the  estimation  of  the 
associated  heights.  Both,  radar  polarimetry  and  radar  interferometry  are  phase  sensitive  techniques. 
The  use  of  polarimetric  SAR  data  has  been  widely  addressed  in  the  last  decade.  The  tight  relation 
between  natural  media  physical  properties  and  their  polarimetric  features  leads  to  highly  descriptive 
results  that  can  be  interpreted  by  analyzing  underlying  scattering  mechanisms.  Interferometric  data  on 
the  other  hand  provide  information  concerning  the  coherence  of  the  scattering  mechanisms  and  can  be 
used  to  retrieve  observed  media  structures  and  complexity.  The  complementary  aspect  of  polarimetric 
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and  interferometric  information  leads  to  a  combination  of  both  approaches.  In  ‘Polarimetric- 
Interferometric  Synthetic  Aperture  Radar  (POL-INSAR)  Imaging’  it  is  possible  to  recover  textural  and 
spatial  properties  simultaneously.  This  includes  the  extraction  of  ‘Digital  Elevation  Maps  (DEM)’ 
from  either  ‘fully  Polarimetric  (scattering  matrix)  or  interferometric  SAR  image  data  takes’  with  the 
additional  benefit  of  obtaining  co-registered  three-dimensional  ‘POL-IN-DEM’  information. 

Over  the  next  few  years  several  free-flying  remote  sensing  satellites  will  be  deployed  in  orbit, 
providing  the  international  scientific,  commercial  and  military  communities  with  a  wealth  of  new  data. 
Many  of  these  will  carry  advanced  multi-channel  imaging  radars  designed  to  combine  various  levels 
of  polarisation  diversity  with  radar  interferometry. 

Polarimetric  interferometry  has  proved  to  be  a  valuable  tool  for  many  applications.  One  is  remote 
sensing  where  it  has  been  shown  in  several  recent  publications  that  by  using  interferograms  in  multiple 
polarisation  channels,  estimation  of  vegetation  height,  underlying  ground  topography  and  mean 
extinction  is  possible  [11,14,23]. 

Polarimetric  Synthetic  Aperture  Radar  interferometry  (POLInSAR)  can  also  be  used  to  enhance  the 
detection  of  military  targets  hidden  beneath  foliage.  The  key  idea  is  to  note  that  for  random  volume 
scattering  the  interferometric  coherence  is  invariant  to  changes  in  wave  polarisation.  On  the  other 
hand,  in  the  presence  of  a  target  the  coherence  changes  with  polarisation.  It  can  be  shown  that  under 
general  symmetry  constraints  this  change  is  linear  in  the  complex  coherence  plane.  These  observations 
can  be  used  to  devise  a  filter  to  suppress  the  returns  from  foliage  clutter  while  maintaining  the  signal 
from  hidden  targets. 


Radar  Polarimetry 


As  shown  in  a  previous  lecture,  an  important  extension  to  single-channel  SAR  remote  sensing  is  the 
utilisation  of  polarised  waves.  A  polarimetric  SAR  system  measures  the  electric  field,  backscattered 
by  the  scene,  including  its  polarisation  state.  The  interaction  of  the  transmitted  wave  with  a  scattering 
object  transforms  its  polarisation. 

One  special  characteristic  of  SAR  polarimetry  is  that  it  allows  a  discrimination  of  different  types  of 
scattering  mechanisms.  This  becomes  possible  because  the  observed  polarimetric  signatures  depend 
strongly  on  the  actual  scattering  process.  In  comparison  to  conventional  single-channel  SAR,  the 
inclusion  of  SAR  polarimetry  consequently  can  lead  to  a  significant  improvement  in  the  quality  of 
data  analysis.  Certain  polarimetric  scattering  models  even  provide  a  direct  physical  interpretation  of 
the  scattering  process,  allowing  an  estimation  of  physical  ground  parameters  like  soil  moisture  and 
surface  roughness  [1 1],  as  well  as  unsupervised  classification  methods  with  automatic  identification  of 
different  scatterer  characteristics  and  target  types  [4,5]. 

SAR  polarimetry  additionally  offers  some  limited  capability  for  separating  multiple  scattering 
mechanisms  occurring  inside  the  same  resolution  cell  and  can  be  deemed  as  a  first  step  in  resolving 
the  ambiguous  scattering  problem  in  SAR,  as  mentioned  above.  With  polarimetric  decomposition 
techniques  a  received  signal  can  be  split  into  a  sum  of  three  scattering  contributions  with  orthogonal 
polarimetric  signatures.  This  can  be  used  for  extracting  the  corresponding  target  types  in  the  image, 
even  in  the  case  that  they  are  occurring  superimposed.  Also,  if  a  signal  is  disturbed  by  undesired 
orthogonal  contributions,  in  this  way  the  relevant  components  can  be  extracted. 
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The  Phase 


In  radar  polarimetry  [1]  we  analyze  the  shape  of  the  transmit  (and  receive)  polarisation  ellipse,  as 
shown  schematically  in  Figure  la,  for  the  purposes  of  improved  information  extraction.  Figure  la 
shows  the  spatial  helix  resulting  from  a  combination  of  horizontal  (H,  in  green)  and  vertical  (V,  in 
blue)  transmitted  components. 


Figure  la  :  The  Polarisation  Ellipse  and  Spatial  helix  decomposed  into  orthogonal  components  x 
(horizontal  H)  and  y  (vertical  V)  [Pottier,  lectures  on  polarimetry] 


By  controlling  the  relative  amplitudes  we  can  rotate  the  polarisation  from  H  through  45  degrees  to  V. 
However,  by  adjusting  the  relative  timing  (phase)  of  the  blue  and  green  components  we  can  also  adjust 
the  shape  of  the  ellipse  as  shown  in  Figure  lb.  It  is  this  combined  amplitude  and  phase  dimension  that 
leads  to  increased  information  content  in  remote  sensing  applications,  since  the  level  of  scattering  we 
observe  from  natural  terrain  depends  on  the  shape  of  this  ellipse.  [1,2, 3, 4, 5] 
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To  represent  this  combined  amplitude  and  phase  control  mathematically,  we  describe  the  wave  using  a 
pair  of  complex  numbers,  ex  and  ey  as  shown  in  equation  1.  The  phase  difference  i.e.  cp  =  arg(exey*) 
then  controls  the  shape  of  the  ellipse,  with  linear  polarisations  defined  by  cp  =  0.  Note  that  the  ellipse  is 
actually  a  dynamic  quantity,  being  the  time  locus  of  the  helix  in  a  fixed  spatial  plane.  Consequently  the 
locus  can  move  clockwise  or  counter-clockwise  (when  viewed  in  the  -z  direction),  corresponding  to 
what  are  termed  left  and  right-handed  polarisations  respectively.  The  set  of  all  possible  left  and  right 
handed  ellipses  can  then  be  conveniently  mapped  onto  the  northern  and  southern  hemispheres  of  the 
Poincare  sphere.  Figure  lc  [1,6] 


Furthermore,  our  ability  to  extract  quantitative  information  rests  on  the  stability  of  this  phase  or  the 
robustness  of  the  spatial  helix  to  small  time  and  spatial  shifts.  A  generic  way  to  define  the  stability  of 
this  helix  is  to  use  the  wave  coherency  matrix  [J],  taken  at  a  point  in  space  with  position  vector  r.  and 
formed  as  an  average  of  all  possible  complex  products  between  ex  and  ey  as  shown  in  equation  1,  [6] 


^  A  ~m — j 

exx  +  e  y  =>  £  = 


e 

Leyj 


[J]  =  (E.E 


(ex  (A*  (-))  (eAK(d) 

_{ey(ryx(rj)  (^(r>;(r))_ 


(1) 
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As  its  name  suggests,  this  matrix  allows  us  to  calculate  not  only  the  wave  intensity  (from  the  diagonal 
components)  but  also  the  coherence,  which  is  a  measure  of  the  phase  stability  of  the  wave,  as  defined  in 
equation  2  [7]. 


rxy  = 


rfcKfcX 


\r„  pi 


(2) 


A  key  benefit  of  employing  ratios  such  as  equation  2  is  that  absolute  amplitude  terms  cancel,  so 
removing  some  of  the  structural  dependence  in  scattering  from  random  media.  It  is  this  observation  that 
shifts  interest  in  polarimetry  towards  the  study  of  ratios  as  potentially  more  robust  indicators  of 
physical  structure  (see  examples  in  table  I). 

In  active  microwave  sensing  we  assume  \^xy  |  =1  for  the  transmitted  wave  and  hence  the  transmitted 

spatial  helix  is  very  stable.  However,  when  the  wave  is  scattered  or  reflected  from  natural  media,  its 
phase  and  amplitude  will  in  general  be  modified  (as  shown  schematically  in  Figure  2a).  This  process 
again  must  be  described  by  a  set  of  complex  numbers,  this  time  by  a  set  of  four,  being  the  elements  of 
the  coherent  scattering  matrix  [S]  defined  as  shown  in  equation  3.  This  matrix  characterises  all  possible 
phase  and  amplitude  changes  due  to  copolar  (diagonal  elements)  and  cross-polar  (off  diagonal) 
scattering.  In  practice,  for  the  common  case  of  backscatter,  the  reciprocity  theorem  for  electromagnetic 
waves  reduces  this  set  to  three  complex  numbers,  as  the  cross-polarisation  terms  are  equal  SHv  =  SVH. 
Note  that  while  this  is  widely  true,  there  are  a  few  special  but  important  cases  where  it  breaks  down,  as 
for  example  in  low  frequency  radio  wave  propagation  through  the  ionosphere,  where  the  earths 
magnetic  field  lines  break  this  reciprocity  symmetry  and  as  a  result  the  cross  polarisation  terms  are  no 
longer  equal.  This  observation  can  be  used  to  calibrate  the  effects  of  Faraday  rotation  due  to  trans- 
ionospheric  propagation,  an  important  issue  for  the  deployment  of  low  frequency  space-borne  radars 
[8]. 


ikr 


£,-—[S]£  =>[*]- 


-V" 

V  7r  - 

'  sHH  " 

V2 SHV 
_  Sw  _ 
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Sw_ 

.  r  iv 
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SHy  =$VH 

[C]  =  )k.k 


(3) 


One  key  idea  in  polarimetry  is  that  if  we  know  all  four  of  these  [S]  matrix  elements  then  we  can 
calculate  the  phase  stability  of  the  scattered  signal  for  arbitrary  incident  ellipse,  using  a  3  x3  covariance 
matrix  [C]  as  shown  in  equation  3.  In  this  way  we  don’t  have  to  actually  change  the  shape  of  the 
transmit  ellipse  (which  would  call  for  control  of  the  antenna  and  microwave  electronics)  but  can 
simulate  the  same  effect  off-line  in  the  processing  stages.  For  this  reason  there  has  been  a  lot  of  interest 
in  the  development  of  microwave  switching  systems  that  are  capable  of  measuring  all  four  elements  of 
[S]  (the  simplest  is  to  switch  each  transmit  pulse  between  X  and  Y  orthogonal  polarisations  with 
simultaneous  reception  of  the  X  and  Y  components).  Note  that  one  important  step  is  to  calibrate  system 
distortion  effects  due  to  crosstalk  (which  causes  problems  with  estimation  of  the  off-diagonal  elements 
of  [S])  and  channel  imbalance  due  to  phase  and  amplitude  distortions  of  the  radar  system  itself.  The 
development  of  robust  calibration  procedures  has  been  a  key  enabling  step  in  the  quantitative 
exploitation  of  this  technology  [9].  Such  systems  are  called  ‘quadpol’  as  they  measure  4  complex 
numbers  for  each  pixel  in  the  image  and  allow  the  user  to  explore  the  whole  Poincare  sphere.  There  are 
currently  several  mature  airborne  quadpol  radar  sensors  with  such  a  capability,  but  significantly  there 
will  soon  be  a  new  generation  of  free-flying  satellite  radars  operating  in  this  mode.  The  European 
Terrasar-X/L,  Japanese  ALOS-PALSAR  and  Canadian  Radarsat-2  are  important  examples.  The  main 
question  then  becomes,  how  can  we  find  the  best  polarisation  combination  to  derive  information 
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products  exploiting  the  scattering  of  waves  from  surfaces  and  vegetation?  To  answer  this  we  must  look 
more  carefully  at  equation  2  and  the  whole  issue  of  coherence. 


Coherence  and  Entropy 


To  calculate  polarimetric  coherence,  we  first  choose  a  pair  of  polarisations  x  and  y,  then  measure  the 
(complex)  components  of  the  signal  in  these  two  channels  and  estimate  the  coherence  by  averaging. 
However,  even  for  a  fixed  wave,  the  coherence  obtained  with  this  method  will  depend  on  the  choice  of 
our  reference  pair  x  and  y  (e.g.  choosing  x=y  will  give  a  coherence  of  1,  while  less  obvious  but  more 
important  is  the  idea  that  for  every  wave  we  can  choose  an  orthogonal  pair  x  and  y  so  that  the 
coherence  is  zero).  This  goes  against  the  idea  that  the  spatial  helix  is  somehow  independent  of  the  co¬ 
ordinates  we  use  to  represent  it,  and  that  consequently  we  should  be  able  to  describe  its  stability  in  co¬ 
ordinate  invariant  terms.  One  way  to  do  this  is  to  describe  the  helix  stability  using  a  generalised 
coherence  or  entropy  (another  popular  way  is  to  use  the  degree  of  polarisation  [6]).  The  wave  entropy 
is  formally  defined  from  the  ratio  of  eigenvalues  of  [J]  (see  equation  4)  and  has  a  value  of  0  when  the 
helix  is  perfectly  stable  and  1  when  it  becomes  noise  like  [2,4,6]. 


0  <H 


w 


(4) 


By  extension,  we  can  also  describe  the  loss  of  helix  stability  after  scattering  by  the  entropy  of  the  3  x3 
covariance  matrix  [C]  in  equation  3,  as  defined  in  equation  5  [2,4] 


(5) 


It  is  important  to  realize  that  this  scattering  entropy  is  characteristic  of  the  scattering  medium  itself.  For 
example,  for  low  frequency  volume  scattering  from  a  cloud  of  ellipsoidal  particles  of  dielectric 
constant  sr  and  axial  ratio  m,  the  normalized  eigenvalues  of  [C]  can  be  evaluated  explicitly  as  shown 
in  equation  6  [2,4,16,17] 


m  =  < 


>  1  prolate  particles 
1  spherical  particles  =>  R  = 
<  1  oblate  particles 


mer  +  2 
m  +  sr+ 1 


\  =2R2  +6R  +  1 
<A2  =  (R- 1)2 

^=(R-\f 


(6) 


For  spherical  particles  (R  =1)  this  leads  to  zero  entropy  but  for  a  cloud  of  ‘wet  dipoles’ (m  and  sr  large) 
the  entropy  rises  to  0.95.  Hence  a  measurement  of  entropy  relates  to  information  about  composition  of 
the  volume.  Importantly,  we  can  estimate  scattering  entropy  numerically  on  a  pixel-by-pixel  basis  from 
quadpol  radar  imaging  data.  Figure  2b  shows  an  example  of  the  entropy  or  phase  stability  of  a  mixed 
scene,  being  the  Oberpfaffenhofen  area  as  collected  by  the  DLR  L-Band  ESAR  system.  We  note  that 
over  non- vegetated  surfaces  (left  bottom  corner  in  figure  2b)  the  entropy  is  low  and  hence  the  scattered 
wave  helix  is  very  stable  for  all  types  of  transmit  polarisation.  This  can  be  exploited  for  quantitative 
moisture  and  roughness  estimation  of  non-vegetated  land  surfaces  by  choosing  appropriate  robust  ratios 
of  scattering  elements  as  shown  for  example  in  table  I.  [10,11,12,13,14,15] 
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The  urban  areas  (upper  right  comer)  in  figure  2  show  moderate  entropy,  but  the  worst  case  arises  for 
vegetation  (lower  right  comer).  Here  we  see  high  entropy  due  to  volume  scattering  by  the  random 
components  of  the  vegetation  cover  (as  in  equation  6).  These  observations  are  independent  of  the  actual 
scene  considered  and  hence  have  been  suggested  by  several  authors  as  suitable  for  robust  unsupervised 
classification  of  land  cover  [3,4,5,13]. 

While  useful  for  classification  and  limited  composition  studies,  such  high  entropy  for  vegetation  cover 
restricts  our  ability  to  fully  exploit  polarisation  for  quantitative  parameter  estimation.  Yet  vegetation 
cover  is  of  prime  importance  in  remote  sensing  applications.  Somehow,  in  order  to  proceed,  we  have  to 
find  a  way  to  reduce  the  entropy.  Importantly  this  can  be  achieved  by  combining  polarimetry  with 
interferometry,  to  form  the  new  topic  of  imaging  polarimetric  interferometry  or  POLInSAR  as  we  now 
show. 


Figure  2a:  Depolarisation  and  Entropy: 


Figure  2b  :  Entropy  image  L-Band  Oberpfaffenhofen,  Germany  DLR  ESAR  Data 
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Controlling  Entropy  :  Volume  Decorrelation  in  Radar  Interferometry 


Radar  interferometry  employs  spatial  separation  by  a  baseline  vector  b  of  multiple  sensors  (for  single¬ 
pass)  or  a  single  sensor  at  multiple  times  (for  repeat-pass)  [18].  It  then  uses  phase  difference  as  a  proxy 
for  elevation,  enabling  determination  of  scatterer  height,  hence  leading  to  products  such  as  high 
resolution  digital  elevation  model  (DEM)  generation.  Again  however,  the  accuracy  of  this  process  is 
governed  by  phase  stability  or  coherence.  In  this  case  we  can  define  a  coherency  matrix  as  shown  in 
equation  7 


hi 


{pMp'M) 

(px(r)p'x(r  +  b)) 


{px(r)p'x(r  + b)} 
(px(r  +  b)p*x(r  +  b)\ 


(7) 


where  ‘x’  corresponds  to  a  single  selected  polarisation  channel.  The  presence  of  vegetation  is  now 
modelled  as  a  finite  bounded  vertical  random  distribution  of  scatterers  with  a  spatial  weighting  to 
account  for  the  fact  that  scatterers  deeper  in  the  volume  will  have  a  smaller  influence  due  to  wave 
extinction.  With  this  model,  the  coherence  of  vegetation  can  be  expressed  as  shown  in  equation  8 
[19,22] 
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where  Bn  is  the  normal  component  of  the  baseline  to  the  line  of  sight.  There  are  two  key  features  of  this 
model: 


•  Coherence  (and  therefore  entropy)  can  now  be  controlled  by  selecting  the  baseline  Bn. 

•  The  interferometric  coherence  is  independent  of  ‘x’  i.e.  of  polarisation 

The  first  means  that,  unlike  in  polarimetry,  we  can  now  design  the  sensor  to  control  the  observed 
entropy  of  vegetation  scattering  (contrast  equations  6  and  8).  However,  the  second  seems  to  indicate 
that  we  do  not  need  polarisation  diversity,  as  equation  8  does  not  change  with  ‘x’.  Why  then  do  we 
need  to  consider  POLInSAR?  The  answer  to  this  apparent  contradiction  is  hidden  in  equation  8  itself. 
We  see  that  the  coherence  is  a  function  of  several  parameters,  the  unknown  height  of  the  vegetation  , 
the  unknown  wave  extinction  and  the  unknown  ground  topographic  phase.  It  follows  that  one  channel 
of  interferometry  by  itself  cannot  be  used  for  unambiguous  parameter  retrieval.  The  situation  is  further 
complicated  by  the  fact  that  for  microwaves  the  extinction  can  be  relatively  small  and  hence  there  can 
be  penetration  of  vegetation  right  down  to  the  underlying  surface.  This  requires  us  to  consider 
combined  surface  and  volume  scattering,  so  forcing  us  to  modify  equation  8  to  at  least  a  two-layer 
model  as  shown  in  equation  9  [19,20,21,22]. 


y  -  JtM  fy  +/4 

•  X 

1  +  74  (9) 

where  /ux  is  the  ratio  of  surface-to-volume  scattering,  which  changes  with  frequency,  vegetation 

density  and  surface  conditions.  However,  it  is  now  that  polarisation  diversity  helps,  as  from  figure  2b 
we  see  that  surface  scattering  has  low  entropy  and  hence  we  can  control  its  influence  in  9  by  changing 
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‘x’  at  the  same  time  as  leaving  the  volume  coherence  unchanged.  Consequently  by  using  POLInSAR 
we  can  increase  the  number  of  observations  faster  than  the  number  of  unknowns  and  hence  achieve 
parameter  estimation  with  a  coherence  or  entropy  under  our  control.  This  is  one  reason  why  there  is 
such  an  interest  in  developing  POLInSAR  sensors  for  vegetation  mapping  [22,23,24].  Several  further 
examples  can  be  found  as  part  of  the  proceedings  of  a  recent  ESA  funded  workshop  focussing  on  this 
topic  (http  ://earth. esa. int/polinsar/) . 

Figure  3  shows  an  example  POLInSAR  product,  obtained  using  the  L-band  airborne  E-SAR  sensor 
operated  by  DLR  in  Germany.  Here  we  show  a  radar-derived  quantitative  tree  height  estimation 
overlaid  on  a  radar-derived  DEM.  It  uses  the  model  of  equation  9  with  polarisation  diversity  over  ‘x’  to 
isolate  the  height  hv  and  cp(z0)  dependence  and  provide  a  map  of  tree  height  over  the  mountainous 
terrain.  Quantitative  comparisons  with  in-situ  measurements  indicate  an  accuracy  of  height  estimation 
around  10%.  [21,22,24].  While  tree  height  is  itself  a  useful  product,  it  can  also  provide  the  basis  for 
various  important  secondary  products.  For  example,  in  Figure  4  we  show  a  forest  biomass  map  derived 
using  the  height  data  in  Figure  3  coupled  to  allometric  equations  derived  from  forestry  tables  for  this 
region  [24].  In  the  upper  Figure  we  also  show  a  conventional  SAR  image  of  the  scene,  which  displays 


Figure  3:  Tree  Height  and  Topography  Estimated  using  L-Band 
DLR  E-SAR  Polarimetric  Interferometric  Data 

none  of  the  important  forest  structural  information  seen  in  the  height/biomass  products.  This  nicely 
illustrates  the  potential  ‘information  gain’  obtained  by  using  POLInSAR  sensors  for  vegetation 
applications. 
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ORGANIZATION 


Figure  4  :  HH  RCS  image  (upper)  and  radar  derived  tree  height/biomass  map 
for  the  same  scene  (lower)  (see  reference  [24]) 


Conclusions  and  Future  Developments 


In  this  treatment  we  have  developed  as  a  theme  the  importance  of  multi-channel  phase  in  radar  remote 
sensing  and  used  it  to  support  the  idea  of  combining  polarisation  diversity  with  interferometry  in  future 
radar  sensors.  Key  to  success  is  the  generalised  coherence  or  entropy  and  key  to  robustness  the 
development  of  physical  models  for  the  interaction  of  polarised  waves  with  natural  surfaces.  We  have 
concentrated  on  one  important  example,  namely  tree  height  and  biomass  estimation,  but  there  are  many 
other  application  areas  where  this  technology  is  being  considered.  Table  I  provides  a  selective  survey 
of  different  geo-physical  parameters  and  examples  of  the  types  of  algorithms  currently  being 
developed.  We  can  see  that  polarimetric  and/or  interferometric  phase  appears  in  every  area.  This  table 
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provides  a  ‘snapshot’  in  time,  each  area  is  ongoing  in  research  and  development  and  exciting  future 
technology  innovations  such  as  bistatic  radar  and  satellite  radar  constellations  will  require  parallel 
improvements  in  our  understanding  of  the  interaction  of  polarised  waves  with  natural  media  in  order  to 
fully  exploit  the  scientific  and  commercial  potential  of  radar  in  remote  sensing. 

By  using  more  than  2  polarimetric  data  sets  polarimetric  interferomerty  can  be  extend  to  more  complex 
approaches.  One  of  possible  approache  is  differential  polarimetric  interferometry.  By  using  3  or  more 
temporal  separated  data  sets  it  could  be  possible  to  enhance  the  potential  of  conventional  differential 
interferometry  by  using  the  polarimetric  information  in  order  to  analyse  the  changes  in  scattering 
processes  over  time.  POL-IN-SAR  imaging,  when  applied  to  ‘Repeat-Pass  Image  Overlay 
Interferometry’,  provides  differential  background  validation  and  measurement,  stress  assessment, 
and  environmental  stress-change  monitoring  capabilities. 


Another  approach  currently  under  investigation  is  polarimetric  SAR  tomography,  which  is  the 
extension  of  conventional  two-dimensional  SAR  imaging  principle  to  three  dimensions.  A  real  three- 
dimensional  imaging  of  a  scene  is  achieved  by  the  formation  of  an  additional  synthetic  aperture  in 
elevation  by  a  coherent  cobination  of  images  acquired  from  several  parallel  flight  tracks.  It  can  be  seen 
as  a  direct  approach  to  resolve  the  SAR  scattering  ambiguity  problems.  The  introduction  of 
tomographic  SAR  offers  the  possibility  of  a  direct  localisation  and  identification  of  all  scattering 
contributions  in  a  volume.  This  greatly  extends  the  potential  of  SAR,  particularly  for  the  analysis  of 
volume  structures  like  for  example  forests  as  shown  in  Figure  5. 


Figure  5  Tomographic  slice  generated  from  13  parallel  flight  tracks  with  a  mutual  distance  of  20m. 
Sensor:  DLR  ESAR  in  L-band.  Scene:  Onberpfaffenhofen,  Germany.  (Reigber[25]) 
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Product 

Radar  Parameter 

Polarimetric  and  Interferometric 
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Table  I :  Examples  of  Geophysical  Parameter  Estimation  using  Radar  Polarimetry  and  Interferometry 
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